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ABSTRACT 


Swell  propagates  across  thousands  of  kilometers  of  ocean  in  almost  unchanged 
parallel  wave  fronts.  Within  the  nearshore  region  however,  refraction  causes  wave  fronts 
to  turn  toward  shallow  depths  transforming  the  wave  field.  The  Nearshore  Canyon 
Experiment  (NCEX)  Pilot,  conducted  from  October  10  to  October  17,  2002,  observed 
wave  transformation  across  the  Scripps  and  La  Jolla  canyon  system  near  San  Diego,  CA. 
Four  Datawell  Directional  Waverider  Buoys,  three  Nortek  Vector  PUV  recorders,  and 
two  pressure  sensors  were  deployed  in  depths  ranging  from  10  to  300  m.  Observed 
energy  density  spectra  and  mean  propagation  directions  were  examined  for  three  case 
studies  representative  of  the  range  of  observed  swell  conditions.  Observations  were 
compared  to  predictions  of  a  back-refraction  model  provided  by  Dr.  William  O'Reilly. 
Observations  indicate  that  refraction  causes  the  waves  to  propagate  along  the  deep  axes 
of  the  Scripps  and  La  Jolla  canyons.  At  the  shallow  canyon  heads,  the  convergence  and 
divergence  of  ray  trajectories  cause  extreme  (2-3  orders  of  magnitude!)  spatial  variations 
in  wave  energy.  Considering  the  complexity  of  the  canyon  environment,  predictions  of 
wave  transformation  agree  surprisingly  well  with  observations. 
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I.  INTRODUCTION 


The  nearshore  environment  is  crucial  to  modem  mval  operations.  The  need  for 
environmental  information  at  a  tactical  level  has  increased  as  naval  power  projection 
continues  to  focus  on  the  littoral  environment.  Amphibious  landings,  nearshore 
reconnaissance,  and  SEAL  team  insertion  are  only  cursory  examples  of  naval  operations 
dependent  on  accurate  nowcasts  and  forecasts  of  the  nearshore  environment. 
Understanding  the  transformation  of  swell  over  complex  coastal  bathymetry;  and  the 
resulting  variations  in  wave  heights,  associated  sediment  transport,  and  surf  zone 
circulation  is  critical  to  the  modern  warfighter. 

Waves  travel  over  thousands  of  kilometers  largely  unchanged  across  ocean 
basins.  Only  a  few  kilometers  from  the  coast,  swell  reach  shallow  water  and  can  be 
dramatically  transformed  by  bottom  topography,  typically  by  processes  such  as 
refraction,  diffraction,  and  reflection.  Refraction,  the  turning  of  wave  trains  caused  by 
the  dependence  of  wave  speed  on  water  depth,  can  be  modeled  precisely  when  high 
resolution  bathymetry  and  boundary  conditions  are  available.  Diffraction  and  reflection, 
which  cause  dispersion  of  wave  energy  away  from  regions  of  discontinuities,  are 
complex  processes  that  are  more  difficult  to  model  accurately  for  realistic  coastal 
settings. 

The  nearshore  environment  off  Scripps  Beach,  near  San  Diego,  California 
provides  an  excellent  location  to  study  the  transformation  of  ocean  swell.  The  relatively 
straight  coastline  extends  along  a  NNW  and  SSE  axis  from  the  Scripps  Institution  of 
Oceanography  to  the  La  Jolla  Beach  and  Tennis  Club  where  the  shoreline  curves  sharply 
to  the  WNW  toward  Pt.  La  Jolla  for  approximately  1000  meters  (Ligure  1).  The  straight 
coastline  is  deceiving,  hiding  complex  undersea  topography  that  causes  extreme  wave 
height  variations  over  distances  less  than  100  m.  The  environment  is  dominated  by  the 
dramatic  depth  changes  of  two  deep  undersea  canyons  that  extend  to  within  250  meters  of 
the  shoreline  (Ligure  1).  The  walls  of  the  larger  and  wider  La  Jolla 
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Figure  1.  Refraction  diagram  for  16- second  waves  from  the  SSW  showing  wave  crests 
and  rays.  Arrowheads  on  the  rays  indicate  the  direction  of  wave  energy  propagation. 
The  shoreline,  and  bottom  contours  are  also  shown,  (from  Munk  and  Traylor,  1947). 


canyon  have  slopes  as  large  as  1:1,  but  even  more  dramatic  is  the  narrow  Scripps  canyon, 
a  northern  branch  of  La  Jolla  Canyon.  The  walls  of  Scripps  canyon  have  been  measured 
to  have  a  slope  as  steep  as  7:1.  Both  rocky  canyons  are  covered  with  a  thin  layer  of  sand 
and  silt,  with  rocky  points  at  the  southern  end  of  the  canyon  rims. 

Munk  and  Traylor  (1947)  conducted  the  first  qualitative  study  of  the  effects  of 

bottom  topography  on  wave  energy  transformation  over  the  Scripps  Canyon.  Refraction 

was  reviewed  in  light  of  other  processes  affecting  the  transformation  of  shallow  water 

waves.  As  waves  approach  the  coastline,  they  begin  to  feel  the  bottom  when  the  water 

depth  is  less  than  one  half  the  wavelength.  Closer  to  the  coast,  when  waves  enter  shallow 

water,  their  wave  speed  equals  J~gli  where  g  is  the  acceleration  of  gravity  and  h  is  the 

water  depth.  Hence,  waves  slow  down  as  the  depth  decreases.  Where  bottom 

topography  is  irregular,  some  parts  of  the  wave  crest  will  feel  the  bottom  before  others. 

Thus,  as  a  wave  crest  approaches  the  shore,  the  portion  of  the  wave  crest  entering  shallow 
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water  slows  down,  causing  the  wave  crest  to  swing  around  toward  areas  of  shallow  water. 
If  the  bottom  slope  is  sufficiently  gentle,  the  wave  crests  tend  to  conform  to  the  bottom 
topography  and  arrive  at  the  coastline  at  near  normal  incidence.  Snell’s  Law  describes 
this  effect:  the  ratio  sin (q)/  c ,  is  conserved  where  q  is  the  angle  of  approach  and  c  is  the 

wave  phase  velocity.  In  areas  with  complex  bottom  topography,  wave  refraction 
diagrams,  such  as  Figure  1,  are  useful  to  visualize  changes  in  wave  direction  and  height. 

In  the  refraction  diagram,  wave  crests  are  shown  as  lines  usually  parallel  to  shore, 
while  the  direction  of  wave  propagation  is  shown  as  a  series  of  rays  that  are  orthogonal  to 
wave  crests.  In  regions  where  rays  converge,  wave  energy  (and  height)  increases,  and  in 
areas  where  ray  traces  diverge,  wave  energy  decreases.  The  example  refraction  diagram. 
Figure  1,  illustrates  the  dramatic  transformation  of  wave  crests  that  are  initially  straight  in 
deep  water.  The  arrows  on  ray  trajectories  indicate  the  propagation  paths  of  incoming 
waves.  The  change  in  spacing  between  rays  indicates  convergence  and  divergence  zones, 
which  in  Figure  1  are  also  signified  by  shading  and  crosshatching.  Notable  features 
include  a  large  area  of  divergence  and  reduced  wave  heights  to  the  south  of  Scripps 
canyon,  extending  to  the  northeastern  side  of  Pt.  La  Jolla.  Convergence  along  the  north 
rim  of  Scripps  Canyon  causes  large  waves  at  Black’s  Beach  on  the  north  side  of  the 
Scripps  Canyon,  and  convergence  along  the  Southern  rim  of  La  Jolla  Canyon  causes 
focusing  at  Pt.  La  Jolla. 

The  wave  period  of  ocean  swell  affects  the  degree  of  wave  refraction.  Waves 
with  longer  periods  "feel"  the  ocean  bottom  in  deeper  water.  Longer  period  swell  (e.g., 
12  and  14  seconds  in  Figures  1  and  2,  respectively)  are  strongly  affected  by  the  canyons 
at  Scripps  Beach,  whereas  shorter  period  wind  waves  do  not  feel  the  bottom  until  they 
have  passed  over  the  canyon  onto  the  relatively  uniform  beach,  and  are  thus  less  affected 
by  refraction  at  this  site. 

Munk  and  Traylor  (1947)  compared  their  theoretical  wave  refraction  estimates 
with  visual  observations  of  wave  heights  collected  by  shore-based  observers  (Figure  2). 
Although  inaccuracies  of  visual  observations,  the  difficulty  of  producing  high  resolution 
refraction  diagrams,  and  the  exclusion  of  the  effects  of  diffraction  or  reflection 
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Figure  2.  Observed  versus  predicted  values  in  wave  height  along  the  beach.  Wave 
heights  were  observed  at  beach  locations  1  through  12.  Predictions  are  based  on 
refraction  diagrams  for  both  14- second  and  12- second  period  swell  from  the  west- 
northwest.  Maximum  wave  heights  are  observed  and  predicted  north  of  the  head  of  La 
Jolla  Canyon,  (from  Munk  and  Traylor,  1947). 


introduced  considerable  uncertainty  in  these  comparisons,  the  observed  and  predicted 
wave  heights  generally  agree  to  within  about  10%.  As  illustrated  in  Figure  2,  refraction 
causes  a  redistribution  of  14  second  swell  wave  energy  at  the  surf  zone  with  variations  as 
high  as  300%  over  only  300  meters  of  coastline.  Actual  variations  may  be  even  higher  in 
the  vicinity  of  the  canyon  heads  where  wave  height  gradients  in  refraction  model  results 
may  have  been  smoothed  out  by  low  resolution  bathymetry  and  where  observations  were 
not  available. 

The  effects  of  swell  refraction  extend  beyond  the  variation  of  wave  height  at  the 
coast.  As  wave  energy  is  redistributed  by  bathymetry,  variations  in  wave  set-up  at  the 
shoreline  cause  horizontal  circulation  cells,  which  in  turn  affect  sediment  transport 
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patterns.  Shepard  and  Inman  (1950)  studied  the  effects  of  wave  refraction  on  nearshore 
water  circulation. 

Observations  using  dye,  drifting  tri- floats,  and  neutrally  buoyant  floats  suggest 
that  the  fundamental  structure  of  nearshore  circulation  results  from  irregularities  in  the 
bottom  topography.  The  observed  nearshore  circulation  cells  can  be  explained 
qualitatively  in  response  to  areas  of  convergent  and  divergent  wave  energy  (Figure  3). 


Figure  3.  Nearshore  circulation  pattern  resulting  from  12.6  s  waves,  (from  Shepard  and 
Inman,  1950). 
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A  major  component  of  nearshore  circulation  cells  are  rip  currents.  Rip  currents 
develop  where  water  flowing  alongshore  away  from  convergence  zones  meets  water 
flowing  from  divergent  zones.  The  combined  water  mass  then  flows  seaward  until  it  is 
deflected  by  larger  longshore  coastal  currents.  Fundamentally,  nearshore  circulation  cells 
depend  on  the  refractive  effects  of  bottom  topography. 

Swell  observed  on  California  beaches  are  the  results  of  storm  events  occurring 
around  the  globe.  In  the  classic  experiment.  Propagation  of  Ocean  Swell  Across  the 
Pacific ,  Snodgrass  et  al.  (1966)  studied  the  evolution  of  swell  from  its  generation  in 
storm  events  in  the  South  Pacific,  through  its  propagation  across  the  Pacific  Ocean,  to  the 
dissipation  of  wave  energy  on  distant  shores.  Analysis  of  their  observations  showed 
surprisingly  little  wave  energy  is  dissipated  in  the  open  ocean,  even  after  swells  have 
traveled  distances  greater  than  12000  nautical  miles.  Figure  4  shows  possible  wave  paths 
that  reach  Southern  California  from  storm  events  in  the  Southern  Hemisphere  and  the 
Gulf  of  Alaska.  In  the  Northern  Hemisphere,  waves  arriving  at  Southern  Californian 
beaches  from  directions  north  of  295°  are  blocked  by  Point  Conception;  and  waves 
originating  from  locations  east  of  150°  are  blocked  by  Point  Eugenia,  Baja  California, 
Mexico.  (Munk  and  collaborators,  1963)  Observations  often  reveal  simultaneous  swell 
arrivals  from  different  propagation  windows. 

The  Southern  California  Bight,  through  which  swell  must  traverse  on  its  path  to 
Pt.  La  Jolla,  contains  many  offshore  islands,  shallow  submerged  banks,  canyons,  and 
other  complicated  coastal  bathymetry,  and  extends  from  32°N  South  of  San  Diego  to 
34.5°N  at  Point  Conception.  Although  islands  shadow  wave  energy  in  their  lee  (Figure 
5),  and  refractive  and  diffractive  effects  decrease  wave  energy,  little  discernable  cross 
swell  is  seen  in  the  lee  of  islands  on  the  spatial  scale  of  the  Southern  California  Bight. 
(Emery,  1958).  Wave  heights,  however,  do  vary  at  different  locations  of  the  Southern 
California  Bight  as  a  result  of  both  the  blocking  effects  of  islands  and  the  scattering 
effects  of  submerged  shoals.  Whereas  refraction  affects  primarily  longer  period  waves, 
island  blocking  affects  all  wave  periods.  Pawka  et  al.  (1984)  compared  observations  of 
wave  height  to  model  results  simulating  the  effects  of  island  blocking  and  wave 
refraction  at  several  locations  across  the  Southern  California  Bight. 
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Figure  4.  Azimuthal  equidistant  projection  centered  on  San  Diego,  California  illustrating 
windows  of  possible  swell  propagation  reaching  the  Southern  California  Bight  from 
storms  around  the  globe.  New  Zealand,  and  Antarctic  pack  ice  significantly  limit  swell 
paths  from  the  storm  centers  in  Southern  Oceans,  (from  Munk  et  al.,  1963). 
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Figure  5.  Schematic  map  of  the  Southern  California  Bight  showing  three  rays  for  16.9  s 
period  waves  incident  to  Torrey  Pines  Beach  near  Point  La  Jolla,  California.  Islands  and 
submerged  shoals  significantly  restrict  possible  swell  paths  arriving  at  Torrey  Pines 
Beach,  (from  Pawka,  et  al.  ,  1984). 


The  model  used  by  Pawka  et  al.  is  a  continuous  spectrum  refraction  model  based 
on  the  conservation  of  wave  energy  and  wave  frequency  along  a  ray.  The  model  requires 
a  linear  wave  field  without  energy  sources  or  sinks,  and  a  slow  variation  of  bathymetry 
and  wave  amplitude.  The  refraction  model  employs  a  continuous  spectrum  to  avoid 
caustic  zones,  and  back-refracts  wave  rays  into  deep  water  spectra.  Rays  intercepting 
islands  are  terminated. 

The  entire  swell  spectrum  is  affected  by  offshore  islands,  sand  banks,  and 
headlands  (Figure  6).  Additionally  diffraction  effects  can  be  important  for  directionally 
narrow  wave  fields  in  the  vicinity  of  complex  bathymetry.  O'Reilly  and  Guza  (1992) 
compared  wave  energy  predictions  from  refraction  and  refraction-diffraction  models  at 
Mission  Beach,  near  San  Diego,  California.  Although  the  refraction  model  is  inaccurate 
when  the  incident  wave  spectrum  is  narrow,  or  when  the  near-field  bathymetry  is 
complex,  in  practical  applications  with  reasonably  broad  spectra  and  mild  bathymetry  the 
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refraction  model  agrees  well  with  the  refraction- diffraction  model.  In  the  case  of 
extremely  broad  directional  incident  spectra,  the  refraction- diffraction  model  produces 
significant  errors  due  to  the  small  angle  approximation  inherent  in  the  model. 
Limitations  of  both  models  in  complex  (e.g.  the  Scripps  Canyon)  environments  where 
both  diffraction  and  large  propagation  angles  are  expected  are  not  well  understood. 
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Figure  6.  Relative  wave  heights  in  the  Southern  California  Bight  predicted  for  October 
16th,  2002  (the  date  of  the  third  case  study)  by  a  refraction- diffraction  model  described  by 
O'Reilly  and  Guza  (1992).  (from  www.cdip.ucsd.edu). 
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D.  EXPERIMENT  AND  ANALYSIS 


Munk  and  Traylor's  1947  pioneering  study  of  wave  refraction  at  Scripps  Beach 
was  primarily  qualitative  owing  to  limitations  of  visual  observations  and  manual  ray 
computations.  Fifty-five  years  later,  computer  systems,  and  advanced  instrumentation 
allow  the  collection  and  processing  of  large  datasets,  and  accurate  modeling  of  wave 
transformation  over  the  nearshore  canyons.  For  8  days,  from  October  10,  2002  through 
October  17,  2002,  data  were  recorded  by  two  pressure  (P)  sensors,  three  Nortek  Vector 
(PUV)  sensors  recording  pressure  and  horizontal  velocity,  and  four  Datawell  Directional 
Waverider  buoys  measuring  wave  height  and  direction.  These  instruments  were 
deployed  in  areas  around  the  Scripps  Canyon  (Figure  7)  as  part  of  a  pilot  study  in 
preparation  for  the  Nearshore  Canyon  Experiment  (NCEX)  scheduled  for  the  fall  of 
2003. 

As  shown  in  Figure  7,  most  instruments  were  deployed  in  the  region  around 
Scripps  Canyon,  with  sparser  coverage  across  the  entire  nearshore  region.  Refraction 
model  computations  were  used  to  select  sensor  locations  in  places  where  the  most 
extreme  wave  energy  transformation  is  expected.  These  extreme  transformations  are 
projected  to  occur  primarily  on  the  north  side  and  at  the  head  of  Scripps  Canyon  where 
the  waverider  buoys  (113-116),  one  pressure  sensor  (P2)  and  one  PUV  sensor  (PUV3)  are 
located.  The  remaining  sensors  (PI,  PUV1,  PUV2)  were  deployed  between  the  canyons 
and  on  the  south  side  of  La  Jolla  Canyon  to  quantify  wave  transformation  over  the  entire 
canyon  region. 

After  a  series  of  bathymetry  survey  sweeps  in  the  vicinity  of  target  pressure 
sensor  locations  verified  acceptable  water  depths  and  bottom  slopes;  two  pressure 
sensors,  designated  PI  and  P2,  were  deployed  on  the  shelf  between  Scripps  and  La  Jolla 
Canyons  (Figure  7).  The  pressure  sensors  were  fitted  in  fiberglass  bottom  tripods  with 
acoustic  releases  for  easy  recovery  (Figure  8).  All  sensors  were  deployed  from  the  R/V 
GORDON  SPROUL  using  differential  GPS  and  targeting  software  for  precise  positioning 
to  within  about  5  meters  of  the  desired  locations. 
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Figure  7.  Locations  of  instruments  deployed  during  the  NCEX  Pilot  experiment  within 
the  nearshore  environment  off  Scripps  Beach.  The  site  features  two  submarine  canyons, 
the  larger  La  Jolla  Canyon  along  a  NNW  to  SSE  axis,  branching  into  the  narrow  Scripps 
Canyon  to  the  east. 
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Figure  8.  Pressure  (P)  sensor  used  in  the  NCEX  Pilot  Experiment.  The  p- sensor  is 
mounted  in  a  red  fiberglass  tripod  equipped  with  lead  feet  to  ensure  stability  on  the  ocean 
bed. 

Divers  from  the  Scripps  Institution  of  Oceanography  measured  the  heights  of  the 
instruments  above  the  sea  bed,  a  measurement  needed  to  accurately  determine  the 
attenuation  of  pressure  over  the  water  column.  Pressure  records  sampled  at  1  Hz  were 
recorded  internally  in  the  instruments.  Recorded  data  were  downloaded  upon  sensor 
recovery. 

Following  the  pressure  sensor  deployment,  and  after  further  bathymetry  surveys, 
the  three  Nortek  Vector  PUV  sensors  were  deployed  (Figure  9).  PUV1  was  deployed 
southwest  of  the  head  of  La  Jolla  Canyon,  PUV2  was  deployed  on  the  shelf  between  the 
canyons,  and  PUV3  was  deployed  near  the  head  of  Scripps  Canyon.  Similar  to  the 
pressure  sensor  moorings,  the  PUV  sensors  were  mounted  in  fiberglass  tripods  to  which  a 
polypropylene  line  and  a  surface  float  were  attached  for  recovery  operations. 
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Figure  9.  The  Nortek  pressure  and  horizontal  velocity  recorder  (PUV)  ready  to  be 
deployed  during  the  NCEX  Pilot  Experiment. 


The  Nortek  Vector  PUV  sensors  contain  a  pressure  sensor  and  a  three- component 
acoustic  Doppler  velocimeter.  The  combined  pressure  (p)  and  horizontal  velocity  (u,v) 
data  provide  height  and  directional  wave  data  equivalent  to  measurements  of  a  pitch  and 
roll  buoy.  The  pressure  and  velocity  time  series  were  recorded  internally  with  a  sample 
frequency  of  2  Hz. 

Waverider  buoys  113  through  116  were  deployed  in  advance  of  the  experiment 
and  continue  to  operate.  These  buoys  and  several  other  permanent  directional  waverider 
buoys  in  the  region  are  maintained  through  the  University  of  California  at  San  Diego's 
Coastal  Data  Information  Program  (CDIP).  The  permanent  Outer  Torrey  Pines  buoy 
designated  Buoy  100,  located  approximately  6.7  miles  west  of  Torrey  Pines  State  Beach 
(Figure  10)  provided  measurements  of  incident  wave  energy  and  direction  spectra  in  deep 
water. 
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The  Datawell  Directional  Waverider  buoys  measure  acceleration  in  three 
directions  and  the  tilt  of  the  buoy  from  which  the  x,  y,  and  z  displacements  are  evaluated 
in  a  fixed  reference  frame.  These  buoys  sample  at  1 .28  Hz,  and  time  series  are  processed 
internally  to  produce  wave  height  spectra  and  directional  moments.  The  buoy  data  are 
transmitted  to  receivers  onshore  via  an  HF  radio  link.  The  buoy  is  moored  with  an  elastic 
line  so  that  wave  motions  at  swell  and  wind  sea  frequencies  (nominally  0.04  -  0.5  Hz)  are 
measured  with  negligible  distortion. 
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Figure  10.  Location  of  NCEX  and  permanently  maintained  directional  waverider  buoys, 
(from  www.cdip.ucsd.edu) 


The  pressure  and  horizontal  velocity  signals  collected  by  the  P  and  PUV-  sensors 
attenuate  with  depth.  Surface  height  spectra  were  estimated  from  the  measured  pressure 
spectra  by  applying  a  correction  in  the  frequency  domain  for  the  hydrodynamic 
attenuation  of  waves  over  the  water  column,  based  upon  linear  wave  theory. 
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Whereas  PI  and  P2  provide  only  surface  height  spectra,  additional  directional 
information  is  available  from  each  PUV  dataset.  Estimates  of  a  mean  wave  direction  at 
each  frequency  were  calculated  using  a  standard  technique  based  upon  pressure- velocity 
cross-spectra  (e.g.  Herbers  et  al.  1999)  Unfortunately  PUV2  did  not  record  meaningful 
data  throughout  most  of  the  experiment,  and  measurements  collected  by  this  instrument 
are  not  included  in  this  analysis. 

Processed  observations  of  the  Datawell  Directional  Waverider  buoys  were 
provided  by  Dr.  William  O'Reilly  of  the  Scripps  Institution  of  Oceanography.  Wave 
height  spectra  were  obtained  directly  from  the  vertical  displacement  measurements 
whereas  the  mean  direction  estimates  at  each  frequency  are  based  on  cross- spectra 
between  vertical  and  horizontal  displacements.  All  estimates  of  wave  spectra  and  mean 
directions  presented  below  are  based  on  3  hour-long  data  records. 
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ID.  OBSERVATIONS 


Observations  made  during  the  NCEX  Pilot  were  compared  in  three  sub -regions 
within  the  nearshore  environment  surrounding  Scripps  Canyon.  These  sub-regions  were 
selected  to  simplify  spectral  comparison  while  facilitating  the  examination  of  wave 
transformation  across  different  paths  of  the  canyon  system.  The  first  sub -region,  north  of 
Scripps  Canyon  includes  the  Waverider  Buoys  113,  114,  and  115,  and  PUV3.  The 
second  sub-region,  surrounding  the  head  of  Scripps  Canyon  includes  buoys  115,  116, 
PUV3  and  P2,  and  the  third  sub-region,  between  Scripps  and  La  Jolla  Canyon,  and  south 
of  La  Jolla  Canyon,  includes  both  pressure  sensors,  PUV1  and  PUV2.  Wave  height  and 
directional  spectra  measured  at  the  offshore  buoy  100  are  included  in  each  sub-regional 
comparison  to  characterize  the  "undisturbed"  incident  wave  field  seaward  of  the  canyons 
and  shelf. 

Three  representative  3-hour  time  intervals  spanning  the  range  of  incident  swell 
conditions  observed  over  the  duration  of  the  experiment  were  selected  as  case  studies:  11 
October  2002  from  1800  to  2100  (PDT),  14  October  2002  from  0300  to  0600  (PDT),  and 
16  October  2002  from  1500  to  1800  (PDT).  These  case  studies  will  be  referred  to  by  the 
date  observations  were  made.  Comparisons  of  energy  and  directional  spectra  are 
restricted  to  the  frequency  range  0.05  -  0.30  Hz,  containing  the  dominant  swell  and  wind 
seas. 

The  incident  wave  field  on  October  14  was  very  similar  to  that  of  October  11,  and 
therefore  this  case  is  not  discussed  here.  On  October  16,  however,  a  significant  shift  in 
the  incident  wave  direction  had  occurred  in  the  swell  band,  and  dramatic  changes  in  wave 
transformations  were  observed. 

A.  11  OCTOBER  CASE  STUDY 

The  incident  wave  conditions  observed  at  buoy  100  (Ligure  11)  show  the 

dominant  swell  with  a  wave  period  of  about  15  seconds  arriving  from  209°  to  214°  True 

North.  This  swell  peak  contains  almost  one  order  of  magnitude  more  energy  than  a 

secondary  peak  observed  at  about  10.5  seconds,  incoming  from  227°.  Higher  frequency 
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local  wind  driven  waves  are  seen  with  periods  less  than  5  seconds,  a  peak  energy  at  4 
seconds  and  a  general  propagation  direction  of  270°.  The  observed  energy  spectra  and 
mean  directions  as  a  function  of  frequency  are  compared  for  each  sub -region  in  Figures 
11  -  13. 

Swell  frequencies  experience  significant  evolution  across  the  canyons.  At  the  15 
second  swell  peak,  energy  levels  are  amplified  (relative  to  the  offshore  buoy  100)  at  both 
PUV3,  located  near  the  most  northern  head  of  Scripps  Canyon,  and  PUV1  at  the  southern 
head  of  La  Jolla  Canyon  (Fig  11-13).  The  focusing  of  energy  at  these  two  locations  is 
consistent  with  the  convergence  of  ray  trajectories  shown  in  Munk  and  Traylor's 
refraction  diagram  (Figure  14).  Before  reaching  the  canyons,  swell  has  been  refracted  to 
the  west  by  the  continental  shelf  south  of  Pt.  La  Jolla.  Additionally,  waves  are  strongly 
refracted  as  they  travel  along  the  axis  of  each  canyon.  An  extreme  example  of  this 
process  is  seen  in  the  directional  spectrum  of  PUV1  (Figure  13),  and  PUV3  (Figure  11- 
12).  Refraction  over  Scripps  Canyon  reverses  the  initial  shift  to  the  west,  whereas  at  La 
Jolla  Canyon,  waves  continue  to  turn  to  more  northern  arrival  angles. 

Swell  energy  decreases  dramatically  between  the  offshore  wave  buoy  100  and  all 
other  nearshore  sensors.  Much  of  this  decay  is  caused  by  refraction  over  the  Pt.  La  Jolla 
shelf.  Additionally,  rays  are  deflected  by  the  two  walls  of  each  canyon,  effectively 
creating  a  'dead  zone'  at  the  head  of  the  canyon.  This  is  evident  in  the  observed  spectrum 
at  buoy  115,  which  records  about  2  orders  of  magnitude  less  energy  than  PUV3,  and  at 
buoy  116  where  the  swell  peak  has  virtually  disappeared  (Figure  12).  The  extreme  decay 
of  energy  at  buoy  115  and  116  can  be  explained  with  the  ray  paths  shown  in  Figure  14. 
Ray  trajectories  diverging  away  from  Scripps  Canyon  carry  very  little  energy,  as  these 
wave  trains  propagate  from  nearly  the  same  offshore  energy  source. 

Buoy  113  measures  energy  levels  closest  to  the  energy  measured  at  buoy  100.  At 
buoy  113's  location,  on  the  shelf  north  of  Scripps  Canyon,  there  is  very  little  convergence 
or  divergence.  The  mean  direction  of  buoy  113  is  more  westerly  than  other  nearshore 
buoys  north  of  Scripps  Canyon.  Wave  trains  refracting  across  the  continental  shelf 
continue  to  refract  westward  as  they  turn  to  shore  normal. 
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Figure  11.  Observed  energy  density  spectra  and  mean  directions  (as  a  function  of 
frequency)  north  of  Scripps  Canyon  on  October  11th.  The  dashed  lines  indicate  the 
incident  wave  conditions  measured  well  offshore  in  deep  water  at  the  Outer  Torrey  Pines 
Buoy,  buoy  100.  Observations  of  wave  conditions  north  of  Scripps  Canyon  are  indicated 
in  color.  Note  the  dramatic  variations  in  wave  conditions  at  swell  (<0.1  Hz)  frequencies. 
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Figure  12.  Observed  energy  density  spectra  and  mean  directions  surrounding  Scripps 
Canyon  on  October  11th.  (same  format  as  Figure  11) 
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Figure  13.  Observed  energy  density  spectra  and  mean  directions  in  the  southern  region 
of  the  canyon  system  on  October  11th.  P- sensor  records  do  not  contain  directional 
information.  Extreme  refraction  is  evident  in  the  observed  turning  of  the  mean  direction 
at  the  swell  peak  by  about  1 10°  from  210°  at  buoy  100  to  320°  at  PUV 1 .  (same  format  as 
Figure  11) 
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P2,  located  between  the  'dead  zone'  east  of  Scripps  Canyon  and  the  convergence 
zone  on  the  shelf  between  the  canyons,  records  about  the  same  level  of  energy  as  buoy 
115.  More  energetic  conditions  in  the  convergence  zone  on  the  shelf  between  the 
canyons  (Figure  14)  are  evident  at  PI  where  the  decay  relative  to  the  incident  wave  field 
is  only  1  order  of  magnitude  (Figure  13). 

Energy  variations  are  smaller  at  the  secondary  peak  with  a  period  of  about  10.5 
seconds.  The  decay  of  energy  levels  relative  to  incident  energy  is  about  a  factor  3-10, 
and  variations  across  the  Scripps  Canyon  are  considerably  smaller  than  those  observed  at 
the  primary  swell  peak. 


Figure  14.  A  refraction  diagram  for  16  second  swell  incident  from  the  SSW  --  conditions 
similar  to  those  observed  during  the  NCEX  Pilot  Experiment,  (from  Munk  and  Traylor, 
1947) 


There  is  very  little  evolution  at  wind  wave  frequencies  (0.15  -  0.3  Hz),  with 
differences  in  energy  levels  less  than  a  factor  2  across  the  canyons,  and  differences  in 
wave  propagation  direction  less  than  10°  across  the  entire  nearshore  region.  The 
relatively  short  wavelengths  of  these  high  frequency  waves  prevents  interaction  with  the 

complex  bathymetry  of  the  canyon  system,  therefore  reducing  energy  and  directional 
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transformations.  At  the  0.25  Hz  wind  sea  peak,  waves  are  transformed  by  refraction  only 
in  depths  less  than  about  10  meters,  well  inshore  of  the  shallowest  pressure  sensor,  PI, 
deployed  in  approximately  16  m  depth. 

B.  16  OCTOBER  CASE  STUDY 

The  October  16  case  study  shows  qualitatively  similar  but  considerably  weaker 
variations  in  energy  levels  (Figure  15-17).  Whereas  swell  energy  levels  measured  on  the 
northern  side  of  Scripps  Canyon  varied  by  about  2  orders  of  magnitude  on  October  11 
(Figure  11)  they  are  within  0.85  orders  of  magnitude  of  each  other  on  October  16  (Figure 
15).  Similarly,  the  amplification  of  energy  levels  at  PUV3  (north  of  Scripps  Canyon)  and 
PUV1  (south  of  La  Jolla  Canyon)  are  less  pronounced  on  October  16,  as  is  the  decay  at 
buoy  116  located  in  the  dead  zone  behind  Scripps  Canyon. 

The  changes  in  energy  variation  are  most  likely  the  result  of  the  shifting  of  the 
incident  swell  direction  to  about  260°.  As  waves  approach  the  canyon  system  from  a 
more  westerly  direction,  they  experience  less  refraction  at  Scripps  Canyon.  Lower 
energy  density  at  PUV3  may  be  due  to  reduced  focusing  at  the  head  of  Scripps  Canyon. 
The  change  in  the  incident  swell  direction  has  exposed  buoy  115  to  more  ocean  swell, 
and  as  a  result  the  wave  field  at  buoy  115  has  a  factor  10  more  energy  than  on  October 
11.  No  longer  in  the  'dead  zone'  caused  by  intense  refraction  at  Scripps  Canyon  on 
October  11,  buoy  116  now  measures  some  swell  energy  that  propagates  across  the 
canyon.  Nevertheless,  buoy  116  measures  about  1  order  of  magnitude  less  energy  than 
buoy  115,  suggesting  that  the  canyon  is  a  strong  barrier  for  these  more  westerly  swell 
arrivals. 

The  energy  amplification  is  also  slightly  weaker  at  PUV1  on  October  16.  With  a 
mean  direction  of  320°  that  suggests  that  waves  continue  to  experience  considerable 
refraction  at  La  Jolla  Canyon.  Similarly  energy  levels  at  PI  and  P2  are  less  reduced 
compared  to  the  open  ocean  swell  energy  (Figure  17). 
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Figure  15.  Energy  and  directional  spectra  at  locations  north  of  Scripps  Canyon  during 
the  October  16th  case  study,  (same  format  as  Figure  11) 
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Figure  16.  Energy  and  directional  spectra  at  locations  surrounding  Scripps  Canyon 
during  the  October  16th  case  study,  (same  format  as  Figure  11) 
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Figure  17.  Energy  and  directional  spectra  at  locations  between  Scripps  Canyon  and  La 
Jolla  Canyon,  and  south  of  La  Jolla  Canyon  during  the  October  16th  case  study,  (same 
format  as  Ligure  11) 
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IV.  COMPARISON  OF  OBSERVED  AND  PREDICTED  SPECTRA 


Accurate  model  predictions  of  wave  energy  and  direction  across  the  nearshore 
environment  are  increasingly  valuable  to  naval  strategists  and  tacticians  planning 
nearshore  naval  operations.  Existing  operational  models  (e.g.  WAM,  WaveWatch, 
SWAN)  can  be  configured  to  predict  the  effects  of  refraction,  but  often  lack  the 
resolution  to  resolve  bathymetry  in  complex  nearshore  environments,  adding 
considerable  uncertainty  to  nearshore  forecasts.  Higher  resolution  models  that  are 
currently  being  assessed  for  use  by  the  U.S.  Navy  require  extensive  verification  with 
experimental  groundtruths.  In  this  chapter,  the  accuracy  of  a  high- resolution  linear  back- 
refraction  model  is  examined  through  comparisons  with  spectra  collected  by  wavebuoys, 
PUVs,  and  pressure  sensors  deployed  in  the  NCEX  Pilot  Experiment  (described  in  the 
second  and  third  chapters). 


A.  THE  BACK-REFRACTION  MODEL 

The  lineal-  back-refraction  model  is  based  upon  the  assumption  that  wave  energy 
is  conserved  along  ray  trajectories,  neglecting  energy  sources  or  sinks.  This 
approximation  works  well  in  small  coastal  regions  where  wave  generation  and  dissipation 
often  can  be  neglected  and  an  accurate  propagation  scheme  is  needed  to  resolve  the 
refraction  of  swell.  Energy  conservation  is  expressed  as 

S(f,Q))=  So(/,  r<y,6))  (1) 

kn  Cg 

where  So  is  the  incident  wave  spectrum  obtained  from  measurements  in  deep  water,  S  the 
predicted  wave  spectrum  at  a  nearshore  location,  /  is  the  frequency,  ko  and  k  are  incident 
and  local  wave  numbers,  &o  and  Cg  are  incident  and  local  group  speeds,  and  V  is  the 
inverse  direction  function  that  can  be  obtained  by  back-refracting  rays  from  the  nearshore 
location  to  the  open  ocean,  (see  O'Reilly  and  Guza,  1991;  and  references  therein)  In 
order  to  accurately  transform  a  spectrum  from  deep  water  to  a  nearshore  instrument 
location,  a  large  number  of  rays  must  be  back- refracted  to  accurately  define  Y(f, 0).  Ray 
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paths  are  initially  computed  using  the  well-known  geometrical  optics  relations  from  a 
selected  nearshore  location  at  1°  increments,  then  additional  ray  paths  are  added  until  rays 
reach  the  deep  ocean  at  1°  increments,  or  the  nearshore  spacing  falls  below  .01°  (O'Reilly 
and  Guza,  1992).  The  predicted  energy  levels  at  each  nearshore  location  are  calculated 
by  applying  Eq.  1  to  the  deep  ocean  incident  spectra  approximated  by  buoy  100. 
Although  many  rays  are  needed  to  adequately  resolve  the  inverse  direction  function,  T, 
over  the  complex  bathymetry,  these  computations  are  performed  only  for  the 
instrumented  sites,  and  thus  the  back-refraction  model  is  computationally  efficient. 
(O'Reilly  and  Guza,  1991) 


B.  MODEL-DATA  COMPARISONS 

Spectra  predicted  by  the  back-refraction  model  are  compared  to  the  observed 
spectra  analyzed  in  Chapter  3.  In  general,  the  model  reproduces  the  observed  energy 
variations  well,  and  considering  the  complexity  of  the  canyon  bathymetry,  the  agreement 
is  surprisingly  good. 

The  model  predictions  accurately  reproduce  the  extreme  variations  of  energy 
density  spectra  and  mean  directions  at  swell  frequencies  observed  on  October  1 1  (Figure 
18).  The  largest  discrepancies  are  seen  at  sites  with  relatively  low  energy  (e.g  at  buoy 
115  and  buoy  116).  At  high  frequencies,  where  waves  are  not  influenced  by  bottom 
topography,  the  predicted  nearshore  spectra  are  virtually  unchanged  from  the  incident 
offshore  spectrum,  in  good  agreement  with  observations. 

On  October  16,  the  predicted  spectra  confirm  that  the  weaker  variations  of  swell 
energy  for  westerly  incidence  angle  (compared  with  the  southerly  swell  arrival  on 
October  11)  are  consistent  with  refraction  predictions  (compare  Figure  18  and  Figure  19). 
However,  a  large  discrepancy  between  the  predicted  and  observed  spectra  is  noted  at 
buoy  115.  Here  the  predicted  swell  energy  is  almost  an  order  of  magnitude  lower  than 
the  observed  energy.  Interestingly,  the  disagreement  is  seen  only  in  the  energy  spectrum, 
whereas  the  observed  and  predicted  mean  directions  agree  well  over  the  energetic  band 
(0.06  -  0.13  Hz)  of  the  spectrum. 
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Figure  18  Observed  (solid)  and  predicted  (dashed)  energy  spectra  and  mean  direction  as 
a  function  of  frequency  on  October  11.  Buoy  100  provided  the  incident  wa\e  spectra  for 
the  back-refraction  model  calculation.  Instrument  locations  are  indicated  with  color: 
Buoy  100  -  black,  buoy  113  -  red,  buoy  114  -  purple,  buoy  115  -  blue,  and  buoy  116  - 
green. 
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Buoy  115  is  moored  in  30  meters  of  water  with  about  60  meters  of  line  near  the 
north  rim  of  Scripps  Canyon.  Energy  levels  at  locations  within  the  resulting  watch-circle 
may  be  significantly  different  from  the  anchor  location  where  the  model  prediction  was 
evaluated.  Buoy  drift  was  considered  as  a  possible  cause  of  the  large  energy  discrepancy 
at  buoy  115.  To  investigate  this  hypothesis,  the  back-refraction  model  was  mn  for 
locations  60  meters  to  the  west  and  60  meters  to  the  north  of  the  buoy  115  anchor 
location.  These  model  predictions  (Figure  20)  indeed  indicate  strong  model  sensitivity  at 
the  steep  wall  of  Scripps  Canyon.  The  predicted  spectrum  60  meters  west  of  the 
presumed  location  of  buoy  115  is  in  excellent  agreement  with  the  energy  spectrum 
measured  by  buoy  115.  However  estimates  of  the  drift  of  buoy  115  obtained  from  the 
GPS  locations  recorded  at  the  buoy  (Figure  21)  show  that  the  buoy  was  confined  to  an 
area  southeast  of  the  anchor  site.  Thus  buoy- drift  does  not  explain  the  modcfdata 
discrepancy. 

A  more  plausible  explanation  of  the  discrepancy  between  observed  and  modeled 
spectra  is  an  inaccuracy  of  the  model  prediction  resulting  from  either  the  breakdown  of 
refraction  theory  over  steep  bottom  slopes  or  the  limited  spatial  resolution  of  the  model. 
The  refraction  model  assumes  a  slowly  varying  smooth  bathymetry,  and  thus  predictions 
may  be  inaccurate  near  canyon  walls.  Model  results  are  generally  more  accurate  for 
locations  with  gentle  bottom  slopes,  such  as  buoy  113.  (Figures  18  and  19)  By 
attempting  to  predict  energy  transformations  on  the  side  of  Scripps  Canyon,  where  the 
gradient  is  as  steep  as  7:1,  we  are  obviously  pushing  the  limits  of  the  model. 
Furthermore,  the  refraction  model  has  a  grid  resolution  of  about  0.0004  degrees  (two  grid 
cells  in  Figure  21)  and  may  not  accurately  resolve  the  large  wave  energy  gradients  at  the 
canyon  wall. 

Overall  the  agreement  between  the  observations  and  refraction  model  predictions 
is  remarkable  considering  the  complexity  and  steepness  of  the  bathymetry,  and  the  results 
of  this  study  suggest  that  a  refraction  model  with  adequate  resolution  can  provide  robust 
forecasts  of  wave  conditions  in  this  region. 
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Figure  19.  Observed  (solid)  and  predicted  (dashed)  spectra  on  October  16.  Buoy  100 
provided  the  incident  spectra  for  the  back-refraction  model  calculation.  Instrument 
locations  are  indicated  with  color:  Buoy  100  -  black,  buoy  113  -  red,  buoy  114  -  purple, 
buoy  115  -  blue,  and  buoy  116  -  green. 
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Figure  20.  The  spectra  observed  by  buoy  100  and  buoy  115,  along  with  the  modeled 
spectra  at  three  locations:  at  buoy  115,  60  m  north  of  buoy  115,  and  60  m  west  of  buoy 
115.  The  modeled  spectra  60  m  north,  and  60  m  west  of  buoy  115  are  in  better 
agreement  with  the  observed  spectrum  at  buoy  115. 
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Figure  21.  GPS  tracked  locations  for  buoy  115  over  the  course  of  the  experiment. 
Although  modeled  spectra  at  locations  60  m  west  and  north  of  the  buoy  (away  from  the 
canyon  wall)  were  closer  to  observations  made  by  buoy  115,  the  buoy  drifted  primarily  to 
the  southeast  (into  the  canyon). 
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V.  SUMMARY  AND  CONCLUSIONS 


Swell  from  distant  storms  propagates  across  thousands  of  kilometers  of  ocean  in 
the  form  of  straight  and  parallel  long- crested  wave  fronts.  As  swell  enters  shallow  water 
and  interacts  with  bottom  topography,  refraction  bends  these  wave  fronts  toward  shallow 
depths, causing  areas  of  energy  convergence  and  divergence.  In  regions  of  complex 
bathymetry,  such  as  offshore  of  Pt.  La  Jolla  near  San  Diego,  CA,  refraction  may  cause 
extreme  transformation  of  swell. 

The  Nearshore  Canyon  Experiment  (NCEX)  Pilot  was  designed  as  a  preliminary 
study  to  quantify  the  effects  of  refraction  caused  by  the  complex  bathymetry  of  the 
Scripps  and  La  Jolla  canyons.  Four  Datawell  Directional  Waverider  Buoys,  three  Nortek 
Vector  PUV  Recorders,  and  two  pressure  sensors  were  deployed  at  locations  of  large 
expected  energy  transformation  from  10  through  17  October  2002.  The  observations 
made  during  the  NCEX  Pilot  show  dramatic  wave  transformation  over  the  submarine 
canyons,  with  refractive  effects  changing  as  a  result  of  incident  wave  direction  and  local 
bathymetry.  These  observations  are  qualitatively  consistent  with  observations  and 
refraction  diagrams  produced  by  Munk  and  Traylor  in  their  pioneering  1947  study  of 
wave  transformation  over  the  same  region. 

Predictions  of  swell  transformation,  using  a  high-resolution  spectral  back- 
refraction  model  were  compared  to  the  NCEX  Pilot  observations.  In  general,  predictions 
and  observations  are  in  good  agreement.  Large  spatial  variations  of  wave  energy  and 
direction  are  well  captured  by  the  model  with  some  discrepancies  in  regions  of  unusually 
low  wave  energy  and  steep  bottom  slopes. 

The  observed  wave  transformation  and  the  skill  of  the  refraction  model  is 
summarized  in  Figure  22,  showing  observed  and  predicted  significant  wave  heights  (over 
the  swell  frequency  range  0.03-0.10  Hz)  at  each  sensor  location  in  the  Pt.  La  Jolla  study 
region  on  October  11.  Relatively  large  significant  wave  heights  are  observed  and 
predicted  at  the  focusing  regions  near  the  head  of  both  canyons.  In  contrast,  the 
significant  wave  height  in  the  'dead  zone'  to  the  east  of  Scripps  Canyon  is  a  factor  10 
smaller  than  the  wave  height  at  the  head  of  Scripps  Canyon.  Additionally,  the  back- 
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refraction  predictions  capture  the  less  pronounced  focusing  of  waves  between  the  two 
canyons.  Considering  the  complexity  of  the  bathymetry  across  the  entire  nearshore 
region  offshore  of  Pt.  La  Jolla,  the  results  of  the  back- refraction  model  are  of  encouraging 
and  confirm  that  refraction  b  the  most  significant  factor  in  wave  transformation  in  this 
region. 


Minutes  West  of  1 17W  Longitude 

Figure  22.  Observed  and  predicted  significant  wave  height  on  October  11.  Red  and  blue 
numbers  correspond  respectively  to  observed  and  predicted  significant  wave  height.  The 
incident  significant  wave  height  observed  at  buoy  100  was  0.38  m.  Model  predictions  for 
PUV1  were  not  available,  thus  only  the  observed  wave  height  is  shown. 
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